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ABSTRACT 


A sequence of solar models has been constructed In order to Investigate 
the sensitivity of the solar radius and luminosity to small changes In the 
ratio a of the mixing length i to the pressure-scale height Hp throughout the 
solar convective envelope. The basic procedure for determining this sensi- 
tivity was to Impose a perturbation In a within the convective envelope and 
then to follow the resulting changes In the solar radius AR and luminosity 
AL for the next 10° yrs. These calculations gave the following results. 

1) A perturbation in a produces Immediate changes In the solar radius and 
luminosity. Initially AL and Aa are related by AL/L ■ 0.30Aa/a. 2) The value 
of the ratio W ■ Alog R/Alog L Is strongly time dependent. Its value_just 
after the perturbation in a is 6.5 x 10”^. 3) The ratio H ■ (Alog L) 

d Alog R/dt Is much less time dependent and Is a more suitable means for 
relating the changes In the solar radius and luminosity. 4) Both of these 
ratios Imply that for any reasonable change In the solar luminosity the cor- 
responding change In the solar radius is negligible. 

I. INTRODUCTION 


During this workshop there has been much discussion about possible 
changes In the solar radius and luminosity over timescales ranging from a year 
or less to a few hundred years. Because of the keen Interest in this topic 
and because of Its obvious relevance to climatic conditions here on earth It 
Is of considerable Importance to determine the sensitivity of the solar radius 
and luminosity to changes In the Interior structure of the Sun. Knowledge of 
this sensitivity together with observational data on any radius and luminosity 
changes would greatly help In understanding the characteristics of the 
physical processes operating within the solar Interior and, as a result, In 
understanding the Influence which these processes might have on the Sun's 
future behavior. In addition. It Is of considerable Importance to determine 
theoretically the relationship between changes In the solar radius and lumin- 
osity resulting from Interior perturbations, since then observational data on 
one of these changes could be used to estimate the size of the other (ref. 1), 
provided, of course, that the physical process causing the perturbations has 
been properly Identified. 

There are otany ways In which the Interior structure of the Sun might be 
perturbed. In approximately the outer 2 per cent of the Sun's mass the 
outward energy flux Is carried largely by convection. Since convection In 
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these layers is a turbulent process, it is entirely plausible that there may 
be random fluctuations in the efficiency of energy transport due, for example, 
to statistical fluctuations in the number of convective cells or to changes in 
the flow pattern* The convective envelope therefore represents one part of 
the Sun where Interior perturbations might be expected naturally to arise* 

Interior perturbations might also occur within the radiative core* Such 
perturbations would alter both the thermal and hydrostatic structure of the 
Sun* The thermal readjustment induced by such perturbations would take place 
over a Kelvin timescale which for the entire Sun is about 10' yrs (ref* 2)* 

The hydrostatic readjustment, however, would take place on a dynamical time- 
scale which is on the order of minutes (ref* 2) and would therefore manifest 
itself almost Instantaneously as a change in the solar radius* 

Studies of other stars provide some evidence that observable changes in 
the radius can result from perturbations within the core* The pulsation 
period of a class of variable stars known as RR Lyrae stars, found both in 
globular clusters and in the field, can be accurately determined by using 
observations spanning several decades* It has been found that the pulsation 
periods of the RR Lyrae stars typically vary at the rate of a few parts in 10^ 
per century* Such changes in the pulsation period can be readily Interpreted 
as changes in the mean stellar radius from one pulsation period to the next* 
The observed rates of period change considerably exceed the values expected 
from the normal evolution of the RR Lyrae stars - a fact that has proved to be 
a long-standing problem* Recent theoretical studies of RR Lyrae models (ref* 
3) have shown that perturbations within the core of these stars can reproduce 
the observed characteristics of the period changes and can thus offer a 
reasonable solution lor this problem. This result suggests that the radiative 
core of the Sun may also be a likely site for the perturbations responsible 
for any changes in the solar radius and luminosity* 

The objective of the present paper is to give the results of one way of 
perturbing the solar interior, namely, by changing the efficiency of energy 
transport by convection throughout the convective envelope* In computing the 
structure of the solar convective envelope it Is necessary to know the value 
of the convective-temperature gradient, i*e*, the actual temperature gradient, 
at each point. The value of this gradient is determined by the requirement 
that the total energy flux carried by both convection and radiation be equal 
to the actual outward energy flux* The convective gradient can range between 
two limiting values, namely, the adiabatic- and the radiative-temperature 
gradients, depending on the degree of convective efficiency* When convection 
is very efficient, the convective gradient approaches the adiabatic gradient* 
This is normally the situation at higher densities and temperatures when the 
thermal energy content of the convective cells is relatively large* At lower 
densities and temperatures, convection can become quite inefficient, and, as a 
result, the convective gradient becomes significantly superadiabat ic and can 
in fact approach the radiative gradient, which is defined to be the tempera- 
ture gradient that would exist if all of the outward energy flux were carried 
by radiation* 



The calculation of the convective gradient as a function of the physical 
conditions at each point In the convective envelope Is generally done accord- 
ing to the prescription of the mixing-length theory. In this theory the tur- 
bulent convective motions which actually cover a tdde range of scale lengths 
are assumed to be represented by convective cells that travel a cb'<racterlstlc 
length I before dissolving Into the surrounding medium. The mixing leng. i t 
Is the main parameter governing the convective efficiency. An Increase In t 
enhances the convective efficiency, thereby lowering the convective gradient. 
Conversely, a decrease In I reduces the convective efficiency, since the con- 
vective cells then cannot transport their excess thermal energy as far before 
dissipation. Ordinarily the value of I at each point is expressed In terms of 
some scale height such as the pressure-scale height H. (*■ dr/dlnP). In this 
paper we will study the consequences of changing the ratio a (> 1/H ) and 
hence the convective efficiency In the solar convective envelope. 9wo points 
should, however, be kept In mind when considering the following results. 

First, there are other ways In which the properties of the solar convection 
could be altered, and hence this paper examines only one type of convective 
perturbation. Secondly, the perturbation results assume that the mixing- 
length theory adequately determines the structure of the solar convective 
envelope at least as far as small perturbations away from the equilibrium 
structure are concerned. 

In the next section we describe first the unperturbed structure of tlu 
solar convective envelope and then the effects which a perturbation In a has 
on this structure at various times following the perturbation. The changes in 
the solar radius and luminosity resulting from a perturbation In a and the 
relationship between these changes are discussed in sections III and IV, 
respectively. We emphasize In section IV the advantages of using the time 
rate of change of the radius perturbation rather than the radius perturbation 
Itself when relating the radius and luminosity perturbations. Finally, a 
summary of the main points Is provided In section V. 

II. SOLAR OONVECTIVE ENVELOPE 

UNPERTURBED STRUCTURE 


In order to examine the unperturbed structure of the solar convective 
envelope, one must first obtain a solar model with the proper luminosity and 
radius at an age of 4.7 x 10^ yrs following the zero-age main-sequence (ZAMS) 
phase. The properties of a solar model are dependent on the assumed composi- 
tion, l.e., the helium abundance Y and the heavy-element abundance Z, and on 
a. For the present calculations Z was taken to be 0.02. The luminosity of a 
solar model Is particularly sensitive to Y, since changes In Y affect the mean 
molecular weight and hence the hydrostatic structure, leading to a change In 
the central temperature. This In turn alters the rate of hydrogen burning due 
to the strong temperature dependence of the nuclear reaction rates. On the 
other hand, a primarily affects the convective envelope and thus the radius. 
Several trial sequences showed that the values 
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Y - 0.2317, a - 1.70 


( 1 ) 


reproduce Che present solar luminosity and radius with an error of about 0.2 
per cent. Accordingly these values were adopted for the model computations, 
and a standard evolutionary sequence was Chen computed from Che ZAMS phase to 
the present Sun. The ZAMS luminosity and radius of this solar sequence were 
0.723 L, and 0.893 R,. 

The unperturbed structure of Che convective envelope In the present Sun 
Is perhaps best Illustrated by Che behavior of the adiabatic-, convective- and 
radiative-temperature gradients (• dlogT/dlogP) , denoted by V , V and 7^, 
respectively. These three gradients are plotted In Figure 1 as functions of 
Che logarithm of the amount of mass between the surface and the given point. 
Here is Che amount of mass within a distance r from the center of the 
Sun. The convective envelope In this model contains 0.016 H^, corresponding 
to log (M^ - M-) ■ -1.796. For a fully Ionised nondegenerate gas with 
negligible radiation pressure V equals 0.40, and we note that 7 approaches 
this value throughout the lnner*part of the convective envelope, ^i.e., for 
log (M^ - Mj.) > ~ -4. Between log (M^ - Mj.) ■ -11 and -4, 7^ Is depressed due 
to Che Ionisation of hydrogen and the first and second Ionisations of helium. 

The difference between 7 and 7^ Is very small for log (M^ - M^) > 7. This 

"adiabatic*' region contains Che bulk of Che mass within the convective 
envelope. In this region energy transport by convection Is very efficient 
with radiation making only a negligible contribution to the outward flux, 
since 7^ » 7 . Because of this high convective efficiency 7 will not be very 
sensitive to a. Just the opposite Is true In the layers near the surface 
(log (M^ - Mj.) < ~ -9), where the convection becomes strongly superadlabatlc. 
As one goes outward through this "superadlabatlc" region, 7 begins to 
exceed 7^ substantially with the maximum superadlabaClclty being reached at 
log (M^ - My) ■ -10.5. The convective efficiency In the superadlabatlc region 
is low due Co the low density and thermal energy content of Che convective 
cells. In Che layers nearest the surface 7 approaches 7 and hence the 
energy transport there la largely by radiation. The structure of the super- 
adlabatlc region trill be strongly dependent on a. The transition between Che 
adiabatic and superadlabatlc regions occurs around log (M^ - My) ■ -8. 

PERTURBED STRUCTURE 

Before discussing Che quantitative results from detailed solar model com- 
putations It Is worthwhile to mention first some further features of the solar 
convective envelope and to consider the physical reasons for the way In which 
the Sun responds Co a change In a. The superadlabatlc region contains little 
mass and has only a small thermal energy content. As a result, the thermal 
timescale of the region Is quite short, on Che order of 1 day (ref. 4), and 
consequently the superadlabatlc region rapidly readjusts to any change In o. 
Within Che adiabatic region 7 i« nearly Independent of a, and therefore 
changes In a that are conflneS to this region will not significantly affect 
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the structure of the convective envelope. However, a change In a throughout 
the entire convective envelope will alter the boundary conditions at the top 
of the adiabatic region, and this In turn will force the adiabatic region to 
undergo both a dynamical and thermal readjustment. The dynamical response of 
the adiabatic region will restore hydrostatic equilibrium on a timescale of 
minutes and will thus be practically Instantaneous. The contraction (or 
expansion) associated with this dynamical readjustment will release (or 
absorb) gravitational potential energy, thereby perturbing the outward energy 
flux Lj. and causing the adiabatic region to depart from thermal equilibrium. 
The timescale for restoring thermal equilibrium Is on the order of l(r yrs 
(ref. 5). The key point to remember Is that the response of L. and hence the 
solar luminosity Is set, not by the thermal timescale of the adiabatic region, 
but by the much shorter thermal timescale of the superadlabatlc region. 
Therefore one would expect a change In a to show up almost immediately as a 
perturbation In the solar luminosity. 


Let us now outline the sequence of events to be expected if, for example, 
a Increases. After about 1 day the superadlabatlc region will have re- 
adjusted both thermally and hydrostatically. As is well-known from stellar 
model computations, an Increase In a leads to a contraction of the adiabatic 
region and hence to the concomitant release of gravitational potential energy, 
resulting In an Increase In the outward flux 1,^ and thus In the solar 
luminosity. The superadlabatlc region will then expand in order to carry the 
additional outward flux, since this is the normal reaction of a region in 
which energy transport by radiation Is Important (ref. 4). Thus one has a 
situation In.whlch the bulk of the convective* envelope contracts on a time- 
scale of ~10^ yrs while the outermost layers initially expand on a timescale 
of days. Observatlonally this would appear as a sudden increase In the solar 
radius followed by a gradual decrease. A similar sequence of events would 
also occur if a decreases except that all of the perturbations would have 
opposite signs. 


In order to verify the above predictions quantitatively, a sequence of 
solar models was constructed in which the time step between models, which is 
normally set by the nuclear timescale of the core, was gradually reduced to 1 
yr. This choice for the minimum time step was made In order to follow the 
rapid changes expected In the solar radius and luminosity tfhile avoiding the 
numerical difficulties sometimes encountered when even shorter time steps are 
used. At this point In the calculations the value of a was Increased by Aa > 
0.01 throughout the convective envelope, and the subsequent evolution of the 
perturbed solar models was followed for about the next 10^ yrs. After the 
change In a the time step was slowly Increased but was always small compared 
with the timescale on which the perturbations were changing. The size of the 
perturbations resulting from this change in a are very small compared with the 
numerical accuracy of typical solar models. For this reason It was essential 
to maintain a high degree of numerical accuracy and especially to minimize the 
importance of numerical noise during the computations. 
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Vftien constructing a stellar model one usually treats the outermost layers 
differently from the interior. In the outermost layers the stellar structure 
equations are integrated inward from the surface to some interior fitting 
point under the assumption of constant L^. This is equivalent to ignoring any 
changes in the gravitational potential energy, i.e., to assuming thermal 
equilibrium. Given several of these integrations, one can then define the 
outer boundary conditions needed for the interior solution. Inside the 
fitting point the stellar structure equations are replaced by difference 
equations which are then solved by an iterative procedure. In the present 
solar models the fitting point was located at log (M^ - M^) > -6. Since 
according to the previous discussion the thermal timescale of the outermost 
10'’° M of the Sun is very short, our implicit assumption of thermul 
equilibrium in these layers should be justified. About '*’5 integration steps 
baaed on a high-order predictor-corrector procedure were used in computing the 
layers above the fitting point. Interior to the fitting point there were 247 
mesh points of which 88 were in the convective envelope. 

There are many sources of numerical noise which can enter into solar 
model computations. For example, stellar structure programs frequently 
contain iterative procedures for determining the density from the equation of 
state, the degree of ionization from the Saha equations and the superadiaba- 
ticity within convective regions. Tight cun/ergence of these iterative pro- 
cedures as tiell as the iterative procedure Involved in the overall convergence 
of the models was required at all times. In addit .on, no changes were per- 
mitted in either the number or distribution of the mesh points. Such changes 
in the mesh points could introduce spurious perturbations by altering the 
truncation error with which the difference equations represent the basic 
differential equations of stellar structure. Special attention must therefore 
be paid to these as well as a number of other sources of numerical noise if 
reliable results are to be obtained. To insure that numerical noise was not 
important in the present calculations, we constructed an additional solar 
sequence in which the perturbation in a was a factor of 10 greater, l.e., 
ta - 0.10. The only difference was the expected scaling in the size of the 
perturbations by a factor of 10. In particular, the ratio of the pertur- 
bations in the solar radius and luminosity changed by less than 2 per cent. 

Let us now consider some of the quantitative results for the readjustment 
of the solar convective envelope after the perturbation Aa • 0.01. Figure 2 
illustrates the difference in the radius A log r between a perturbed model and 
the basic unperturbed model as a function of M. - within the convective 
envelope. The four curves labelled a. b, c and d correspond to four perturbed 
models having ages of 1, 4900, 49,000 and 310,000 yrs, respectively, following 
the perturbation in a. The contraction is not yet noticeable in oiodel a, 
because the time elapsed since the perturbation has been too short. Moreover, 
the increase in the surface radius A log R in model a due to the expansion of 
the superadlabatic region amounted to only 5 x 10**^. By model d the rate of 
contraction has slowed substantially so that this model is approaching the 
equilibrium structure for the new value of a. The amount of the contraction 
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Is considerably greater nearer the surface, and thus the convective envelope 
does not contract uniformly* This contraction Increases the weight of the 
convective envelope on the radiative core, thereby causing the core also to 
contract, as Indicated In Figure 2 for > 0*016 M^* 

The rate of release of gravitational potential energy e In ergs/gm/sec 
within the convective envelope Is shown In Figure 3 for each*of the four per- 
turbed models plotted In Figure 2* The maximum rate of contraction of the 
convective envelope occurs Immediately after the perturbation In a, and thus 
the largest values of t are produced at this time* However, the radiative 
core does not begin to lontract until after there has been a decrease In the 
radius of the convective envelope and hence a change In the boundary con- 
ditions at the edge of the core* This explains why the release of gravita- 
tional energy In the core Is negligible In model a while It becomes important 
In the later mod'>ls* We note that c Is negative for ~ Mj, < 6 x 10”^ In 
model a due to the expansion of the outer layers of the convective envelope* 

In the present calculations this expansion disappears 1 year after the pertur- 
bation In a; It might actually disappear sooner If shorter time steps are 
used* The slowlng-down of the contraction with time, as Indicated by the 
decrease Is apparent In going from models a to d* 

The release of gravitational potential energy perturbs the outward flux 
Lp at each point within the convective envelope. This flux perturbation 
Is Illustrated in Figure 4, where the difference In between each of the 
four perturbed models In Figure 2 and the unperturbed model Is shown over the 
same Interval In as In Figure 3* The behavior of the flux perturbn- 

tlon In time Is somewhat complicated In the Inner half of the convective 
envelope due to two competing effects. Between models a and c the c^ttribu- 
tlon to the flux psrturbatlon from the contraction of the core lncr^>^«ea, 
while at the same time the contribution of the convective envelope decreaaes. 
The drop In AL^, for - M_ < 6 x 10“^ in model a Is again associated with 
the Initial expansion of the outermost layers* 

The above discussion has focused on the structural readjustment that 
takes place within the solar convective envelope following a perturbation 
In o* We now turn our attention to the question of tdiat potentially observable 
changes a perturbation In a might produce In the solar radius and luminosity* 

III. CHANGES IN THE S0L>R RADIUS AND LUMINOSITY 


The changes in the solar radius Alog R and luminosity Alog L during the 
first 8 X 10^ yrs after the perturbation Ir o are presented in Figure 5. The 
zero-point of the time scale In Figure 5 as well as In all subsequent figures 
corresponds to the time t when the perturbation Aa * 0.01 was imposed within 
the convective envelope* The response of the solar luminosity to this pertur- 
bation appears to be nearly Instantaneous for the time resolution of this 
figure. Following the large InitlaJ response Alog L decays with an e-foldlng 

yrs. By the latest times shown In Figure 


time on the order of a few times 10 
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5 the perturbed solar models ere e^proeching their new equllibrius structure 
which is characterized by a decrease in log R and an Increase in log L* The 
present results demonstrate that the initial response of the solar luminosity 
considerably exceeds the difference In lo^ L between the unperturbed and new 
equilibrium states. Also plotted in Figure 5 is the change in the rate of 
hydrogen burning Alog L„. The contraction of the core, as indicated pre-> 
viously, raises the temperature in the layers near the center, thus increasing 
the rate of the nuclear reactions. 

It is of some interest to examine the behavior of Alog L and Alog R 
immediately following the perturbation in a. Figure 6 shows this behavior for 
Alog L. The time scale in this figure has been expanded by approximately a 
factor of 2000 compared with Figure 5 and consequently covers only the first 
400 yrs after the perturbation in a. Even on this expanded timescrle there is 
a sudden response of the solar luminosity at t • 0. This response would 
actually have been more abrupt if tloK steps less than 1 year had been used in 
the computations. This result confirms our previous conjecture that changes 
in convective efficiency of the type considered here will almost immediately 
affect the surface luminosity. Ue note from Figure 6 that Alog L is nearly 
constant over a timescale of several hundred years. From these results it 
follows that the change in the solar luminosity produced by a perturbation Aa 
is given by 


~ - 0.30 (2) 

L a 

for short times after the perturbation. A similar expression has been derived 
by Dearborn and Blake (ref. 4), who found a coefficient of 0.f4 on the right- 
hanJ side of equation (2). 

The more complicated behavior of Alog R is Illustrated in Figure 7 for 
the same time Interval as in Figure 6. The ordinate in Figure 7 has been 
expanded by roughly a factor of 1000 in comparison with Figure 5. The sudden 
increase of the solar radius due to the expansion of the supersdlabatlc region 
if readily apparent at t « 0. This initial expansion is folljired by an 
overall contraction of the convective envelope and hence in the solar radius 
as the adiabatic region reacts to the change In a. At t •• 250 yrs the radius 
again equals its unperturbed value. The maximum value of Alog R just after 
the perturbation in a was quite small, only 5 x lO"', which explains why the 
initial expansion was not evident in Figure 5. This maximum value of 
Alog R Is related to the perturbation Ao by the equation 

^ - 2.C X 10"^ ^ . (3) 

Two features of Figure 7 should be emphasized. First, the value of Alog R is 
strongly time dependent even over the short time interval covered by this 
figure. Second, the rate of change of Alog R, d Alog R/dt, is, in contrast, 
nearly constant. 
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In this section we have discussed each of the changes Alog L and 
A log R separately. We now wish to consider how the»e changes are related to 
each other. 

IV. RELATIONSHIP BETWEEN CHANGES IN THE SOLAR RADIUS AND LUMINOSITY 


One of the objectives of previous studies (refs. 1, 4, 6, 7, 8) on the 
effects of perturbations In a was to determine the ratio 


W 


Alog R 
A log L 


(4) 


This ratio can be straightforwardly obtained from the present calculations to 
give the results shown In Figure 8, where the time Interval Is the same as In 
Figures 6 and 7. The strong time dependence of W Is Immediately evident. In 
fact, the value of W changes sign at t >■ 250 yrs. Since Alog L Is nearly 
constant over the time Interval In Figure 8, this time dependence Is actually 
a reflection of the strong time variation of Alog R. The values of W In Figure 
8 can be approximated by the equation 


W(t) - 6.5 X lO"^ - 2. : X lO’^ t, 


(5) 


where t Is in years. The original estimates of V (t * 0) ranged from 0.075 to 
5 X 10”^ (refs. 1, 4). More recent determinations have averaged from 5 x 10“^ 
to 10 X 10“^ (refs. 6, 7, 8) and are therefore In agreement with the present 
value . 


One would like to use W to determine, for example, the change In the 
solar luminosity associated with observational estimates for changes In the 
solar radius. However, there are two major disadvantages with using W for 
this purpose. First, It Is only appropriate to use W If the perturbation in 
a has occurred during the time Interval spanned by the radius observations. 
Otherwls . any observed change In log R would actually be the change between 
two pertu..bed states rather than between the unperturbed and perturbed states. 
From the last section we know that a perturbation In a gives rise to changes 
Alog R and Alog L that persist for several times 10^ yrs. Thus, If an 
observed change In the solar radius Is ascribed to a perturbation In a, the 
probability that this perturbation occurred during the Interval of the obser- 
vations Is very small. Second, there Is the problem caused by the strong time 
dependence of W. Even If the Hrst disadvantage is Ignored, one must still 
know how much time has elapsed since the perturbation In a In order to compute 
the proper value of W from equation (5). 

The above difficulties can be overcome by using an alternative expression 
relating Alog R and Alog L, namely, the ratio 

H- L_lAi2SJ: (6) 

Alog L dt Alog L dt 
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This ratio Is plotted In Figure 9, where the tlute Interval Is again the same 
as In Figures 6 and 7* The average value of H In Figure 9 Is 

H - -2.6 X lO’^ yr”^ (7) 


The variation tn the value of H Is substantially less than was the case for W. 
This result lu not surprising In view of our previous comments that Alog L and 
d Alog R/dt In Figures 6 7 are nearly constant. The fact that H Is not 

strongly time dependent can also be given a straightforward physical explana- 
tion. The rate of contraction of the convective envelope d Alog R/dt deter- 
mines the rate of release of gravitational potential energy which In turn 
determines the luminosity change Alog L. Thus d Alog R/dt and Alog L actually 
represent different ways of measuring the same quantity, namely, the mean 
value of e In the convective envelope, and consequently we would expect this 
ratio to b§ approximately constant at least for short times following the 
perturbation In a. Over much longer time Intervals, however, the value of H 
will change significantly, ae Is Illustrated by Figure 10, but even here the 
relatlx'e change Is much less than that shown by W In Figure 8. For example, 
after 10^ yrs the value of H differs by only a factor of 2 from Its value at 
t - 0. When using H to relate Alog L to an observed radius change, one Is 
Implicitly assuming that the perturbation In a occurred prior to the time of 
the observations, but, as mentioned before, this Is very likely to be the 
case. We conclude therefore that the inherent disadvantages of the ratio W 
can be circumvented to a large extent by using the ratio H. 


Dunham et al. (ref. 9) have reported a decrease In the solar radius of 
0.70 ± 0.12 between 1923 and 1980 from measurements of the size of the path of 
totality during a number of solar eclipses. The corresponding change In log R 
is thus -3 X 10“^. Let us now see what this observational result Implies for 
the change In the solar luminosity under the assumption that a perturbation in 
a Is responsible for the radius change. There are two cases to consider. 
First, let us assume that the perturbation In a occurred sometime after 1925 
so that W is the appropriate ratio to use. From equation (5) It follows that 
5.1 X 10“^ < W < 6.5 X lO”"*. The change In log L determined from these values 
of W lies In the range -0.62 < Alog L < -0.49, Implying that the solar 
luminosity In 1925 differed from the present luminosity by a factor of 3 or 
4. As the second case, let us assume that the perturbation In a occurred 
before 1925 so that we must apply the ratio H. The radius measurements then 
give -6 X 10“® yr”^ for the average value of d Alog R/dt since 1925. By 
combining this observational result with the value of H from equation (7), we 
find that Alog L ■ 2.2, again Implying an Impossibly large change In t)<e .olar 
luminosity. The change In log L would have been even greater If a si>.al. 
value of 1H| had been used, as would be appropriate for later times 
according to Figure 10. We conclude therefore that the change in the t^olar 
radius since 1925 either has not been as large as reported by Dunham et al. or 
has been produced by some process other than the one studied in this paper. 
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V 


SUMMARY 


From the present results It is possible to draw the following 
conclusions: 

1) Changes In the efficiency of convection throughout the solar convec- 
tive envelope lead to sudden changes In both the solar radius and luminosity* 
The relationship between the change in the luminosity and the change In a is 
given by equation (2). 

2} The value of the ratio W > Alog R/Alog L Is strongly time dependent. 
For this and other reasons W does not seem to be a very suitable means for 
relating changes In the solar radius and luminosity. Immediately after a 
perturbation In a the value of W Is 6.5 x 10"^. 

3) A more satisfactory way to relate the radius and Itunlnoslty changes Is 
represented by the ratio H ■ (Alog L) ^4 Alog R/dt. This ratio la much less 
tlTO dependent, varying from -2.6 x 10”** to -1.3 x 10“® yr”^ during the first 
10^ yrs following a perturbation In a. 

4) According to the present values of W and H, any observatlonally detect- 
able change In the solar radius would Imply an Impossibly large change In the 
solar luminosity. Consequently changes In convective efficiency of the type 
considered here cannot be responsible for any observed radius changes In the 
Sun. 
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Figure !• Variation of the adiabatic-, convective- and radiative- 

temperature gradients V^, and respectively, as functions 
of log (M^ - M^) within the solar convective envelope 
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Figure 2« Perturbation In the radius Alog r as a function of - Mj. 

within the solar convective envelope at four times following the 
perturbation Aa > 0.01 
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Figure 3* 


Rate of release of gravitational potential energy e as a 
function of within the solar convective envelope at four 

times following the percurhation Aa « 0*01 
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Flgijre km Perturbation In the outward flux AL^ as a function of 

within the solar convective envelope at four times following the 
perturbation Aa - 0*01 
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Figure 6* 
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Time dependence of the change in the luminosity Alog L shortly 
after the perturbation Aa - 0.01 vithin the solar convective 
envelope 




Figure 7. Time dependence of the change In the surface radius Alog R 

shortly after the perturbation Aa - 0.01 within the solar con- 
vective envelope 
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Figure 8* 


TIm dependence of the ratio W - Alog R/Alog L shortly after the 
perturbation Aa * 0«01 vithin the solar convective envelope 
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Figure 9 . Time depeoder.ee of the ratio H ■ (Alog L)”^d Alog R/dt shortly 
after the perturbation Aa - 0.01 within the solar convective 
envelope 
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Figure 10. Tlae dependence of the retio H “ (Alog L) d Alog R/dt 

following the perturbation Aa - 0.01 within the solar convective 
envelope 



